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Abstract: Liquid scintillating phantoms have been proposed as a means to perform real-time 
3D dosimetry for proton therapy treatment plan verification. We have studied what effect radiation 
damage to the scintillator will have upon this application. We have performed measurements of the 
degradation of the light yield and optical attenuation length of liquid scintillator and water-based 
liquid scintillator after irradiation by 201 MeV proton beams that deposited doses of approximately 
52 Gy, 300 Gy, and 800 Gy in the scintillator. Liquid scintillator and water-based liquid scintillator 
(composed of 5% scintillating phase) exhibit light yield reductions of 1.74 ±0.55% and 1.31 ± 
0.59% after « 800 Gy of proton dose, respectively. Whilst some increased optical attenuation 
was observed in the irradiated samples, the measured reduction to the light yield is also due to 
damage to the scintillation light production. Based on our results and conservative estimates of the 
expected dose in a clinical context, a scintillating phantom used for proton therapy treatment plan 
verification would exhibit a systematic light yield reduction of approximately 0.1% after a year of 
operation. 
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1. Introduction and Purpose 

Dose verification Quality Assurance (QA) measurements are critical in proton therapy for the val¬ 
idation of patient treatment plans []II]. These measurements are especially important in Intensity 
Modulated Proton Therapy (IMPT), which scans a focused proton beam spot over the target region 
many times at different energies and intensities to improve dose conformity. IMPT carries in¬ 
creased treatment planning complexity relative to passively scattered treatment modalities, making 
the need for dose verification over all beam energy layers especially critical. 

An emerging concept for proton therapy dose verification systems is the use of liquid scin¬ 
tillator phantoms with position sensitive photon sensors to image the 3D dose distribution with 
the phantom volume. These novel detectors present several advantages over the current QA mea¬ 
surement systems that rely upon the physical movement of a two-dimensional array of gaseous 
detectors within a water phantom [0,0. In particular, liquid scintillator-based systems permit sub¬ 
millimeter spatial resolution 0], and have the potential to perform validation measurements over 
all energy layers in no more time than is required for the patient’s exposure to the beam. This 
real-time dose verification is an especially attractive feature of the active phantom approach, as it 
facilitates a high patient throughput for clinical operations. 

Water-based Liquid Scintillator (WbLS) is a newly-developed water-scintillator emulsion 10 
that is a promising candidate material for next-generation particle physics detectors [0. As WbLS 
is predominantly water, it offers improved water-equivalence for a scintillator-filled dose-verification 
phantom, thereby offering superior dosimetry performance relative to traditional liquid scintillator. 
WbLS also largely mitigates the flammability and toxicity hazards that are traditionally associ¬ 
ated with liquid scintillators. We are currently exploring the possibility of using WbLS for dose 
verification QA in proton therapy. 
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Figure 1. Details of the proton beam irradiation geometry; (a) a schematic drawing of the sample holders 
used in the irradiations, (b) the proton beam intensity profile incident upon the scintillator volumes (outlined 
areas), taken using a phosphor screen coupled to a CCD. The fiducial marks in (b) define the 10 cm by 10 
cm region about the beam centre. 

One aspect of the use of scintillating materials for dose-verification systems that has yet to 
be comprehensively investigated is the potential for radiation damage to affect the scintillating 
performance of the phantom volume. Such damage could lead to systematic measurement errors 
due to a variable scintillation light yield as a function of dose deposited. In this study, we have 
quantified fhe relafionship befween accumulafed dose and scinfillafing lighf yield for pure liquid 
scinfillafor and WbLS samples. Our measuremenfs of changes fo fhe lighf yield have been assessed 
in ferms of fhe expecfed effecfs upon QA defector operation. 

2. Experiment 

2.1 Materials 

For fhe pure liquid scinfillafor samples we used a Linear Alkyl Benzene (LAB) solvenf wifh 3 g/L 
2,5-Diphenyloxazole (PPO) as primary fluorophore and 15 mg/L L4-Bis(2-mefhylslyryl)benzene 
as a wavelengfh shifting fluorophore. Our WbLS was composed of 5% LAB by mass wifh 2 g/L 
PPO in fhe organic phase as fhe primary fluorophore. 

2.2 Irradiation A 

Profon irradiafions were carried ouf af fhe NASA Space Radiafion Laboratory (NSRL) af Brookhaven 
National Laboratory. Duplicate samples in polypropylene vials were mounfed in fwo HDPE mounfs 
and exposed fo a 10 cm x 20 cm beam of 201 MeV protons (figure |T^. The uniformify of fhe beam 
infensify across fhe fargef area was measured fo be approximafely ±13% (figure [T^. Two consec- 
ufive beam exposures were used fo deliver approximafely 300 Gy and a furfher 500 Gy fo samples 
of pure liquid scinfillafor and WbLS. Samples from fhe rear mounf were removed following fhe 
firsf exposure, and samples in fhe fronf mounf were irradiafed by bofh exposures. The dose accu¬ 
mulafed in fhe rear mounf was verified using a Fricke dosimefry solufion fhaf was placed in fhe 
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sample holder along with the scintillating samples. The Fricke measurements suggest that the first 
exposure delivered a dose of 328 ± 10 Gy to the Fricke vial. The 800 Gy accumulated by the front 
holder’s samples is twice the maximum dose measurable with Fricke dosimetry []7|], so the NSRL 
beam instrumentation - a calibrated ionisation chamber placed upstream from the samples - was 
relied upon for the dose estimate. 

During the irradiations, online measurements of the scintillation response of liquid scintillator 
and WbLS samples were taken to search for short term radiation damage effects. Liquid Scintillator 
and WbLS samples were placed 50 cm behind the HDPE holders and measurements were carried 
out using photomultipliers and a digital oscilloscope. A gas ionisation chamber was used as a 
reference to correct for variations in the beam intensity. No short term damage to the light yield of 
either sample was observed, to a precision of 3.3% and 2.4% for the WbLS and liquid scintillator 
samples, respectively. The dose delivered to the samples is estimated to be 556 Gy, using the 
simulation model described below. 

Empty vials were also placed in the beam area and later filled with scintillator and water-based 
liquid scintillator to check for radiation-induced leaching from these vials. Measurements of the 
scintillating materials that resided in these vials for approximately 1 month revealed no change in 
the light yield, suggesting that vial leaching does not affect our results. 


Table 1. The simulated mean energy deposit per interacting proton, relative beam intensity, and relative dose 
in the irradiated samples. All results refer to a vial placed at the beam centre. 


Sample Eocation 

Mean Energy Deposit 

Relative beam intensity 

Relative dose 

Eront HDPE holder 

9.37 MeV 

1 

1 

Rear HDPE holder 

9.75 MeV 

0.957 

0.996 

Online measurement 

9.99 MeV 

0.652 

0.695 


2.2.1 Simulation Model 

We have modeled our irradiation experiments using a Geant4 vlO.O [H] simulation of the measure¬ 
ments. The calculated dose information for the irradiated samples is presented in table [I]. The doses 
delivered to the front and rear vials were identical to within 0.5%. Eigure @ plots the original beam 
position of protons that interact in the rear vial. The relatively sharp image of the liquid volume 
suggests the beam divergence upon entry into the rear vial is on the order of 1 mm. Since the beam 
spreading was small compared to the length scale of the variations in the incident beam intensity 
(figure |I^, the incident beam profile was used to estimate the relative dose at the sample locations 
for the rear sample holder. We attribute the excess events above the scintillating volume in figure ^ 
to delta rays that travel through the air above the scintillator. 

2.3 Irradiation B 

A second irradiation was carried out to investigate the effect of low doses upon the scintillator light 
yield. A beam of 201 MeV protons was used to deliver a dose of approximately 20 Gy and 50 Gy 
in two irradiations to an irradiation geometry similar to figure |T^ A 21.5 cm thick HDPE block 
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Figure 2. The initial beam location of simulated protons that deposit energy in the scintillator encapsulated 
by the central rear vial. The shaded region denotes the location of the scintillator. 

was used to slow the incident proton beam, allowing the vials in the rear holder to receive a greater 
dose from stopping protons. Fricke dosimetry measurements of the front and rear vials suggest 
that the dose ratio in the rear relative to the front was 1.6. A beam profile similar to figure |T^ was 
obtained. 

2.4 Light yield measurements 

Scintillation measurements were carried out using a Beckman LS6500 liquid scintillation spec¬ 
trometer. A *^^Cs source was placed inside the instrument’s lead shielding near the sample cham¬ 
ber to excite the scintillator. The Compton electron spectrum due to this source was measured 
repeatedly for all samples both prior to and after the irradiation. All irradiated scintillating samples 
were dispensed into unirradiated vials and compared with an unirradiated control. Some instru¬ 
mental artifacts were observed in the measured scintillation spectrum. These fluctuations in the 
histogrammed pulse height data exhibited a bin-to-bin correlation in a quasi-periodic fashion (fig¬ 
ure^. The amplitude of these systematic variations did not vary strongly with the number of events 
in the histogram, so that acquiring for longer periods could minimise this effect. The precise cause 
of these artifacts is not known. The effect of the instrumental artifacts is taken into account in the 
assigned systematic uncertainty. 

The change in the scintillation light yield upon irradiation was assessed relative to the control 
sample using an algorithm that compared the sample’s Compton electron spectrum with that of the 
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Figure 3. An example of the measured Compton electron spectrum for the control (black) and sample of 
liquid scintillator irradiated by wSOO Gy of 201 MeV protons (red). The green trace shows the irradiated 
spectrum scaled by the amount of light yield loss as determined by the algorithm described in the text. The 
non-statistical fluctuations are an artifact of the liquid scintillation spectrometer readout. 


control. A multiplicative factor was used to scale the sample spectrum along the channel axis. The 
histogram binning was correctly accounted for in the scaling operation. The change in the mea¬ 
sured light yield was determined as the scaling factor that minimised the between the measured 
and control spectra in a region above a channel threshold and below the bin in the control spec¬ 
trum where the number of counts fell below 20. The optimal channel threshold was determined 
empirically for each sample. The operation of this algorithm is illustrated in figure 0. 

To account for the possibility that our algorithm’s multiplicative factor is a biased estimator of 
the light yield, a Geant4 vlO.O simulation was developed to synthesize a similar Compton electron 
spectrum to that measured in our detector. An isotropic source of 662 keV gamma rays was placed 
2.5 cm from a scintillating volume identical in size to that used in our measurements. A 2.5 cm 
thick cylindrical lead shield with an inner diameter of 7.5 cm surrounded the scintillating volume 
and source. The simulated spectrum of energy deposits in the scintillator was convolved with a 
Gaussian of sufficient width to give an empirical match for the measured liquid scintillator and 
WbLS spectra. Systematic errors in our light yield algorithm were investigated by comparing the 
algorithm’s estimate of the scaling factor with the known scaling factor. This analysis shows that 
the algorithm underestimates the amount of light yield loss for both the pure liquid scintillator 
and the WbLS, with the lower amplitude WbLS spectrum having a larger error. The magnitude 
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(a) Liquid Scintillator 


(b) Water-based Liquid Scintillator 


Figure 4. The optical absorbance spectra for unirradiated samples (black), and samples irradiated with 
protons (blue and red). The average doses received by samples in the same irradiation holder are given. 


of this error was also found to increase as the reduction in the light yield increased. The required 
correction was less than 0.5% of the light yield for all measured light yield losses, and we have 
corrected our reported results for this effect. 

2.5 Optical Characterisation 

Optical absorbance measurements of the irradiated samples were carried out using a Shimadzu 
UV1800 spectrophotometer. We observed an increase in the optical attenuation coefficient of the 
irradiated scintillators that increased with increasing dose delivered to the scintillator (figure 


3. Results 

The reduction of the light yield at the proton doses received by the samples is given in figure ^ 
The samples were measured approximately 24 hours following irradiation to allow short lived acti¬ 
vation products to decay. No increase in the count rate was observed relative to the pre-irradiation 
measurements, suggesting that long lived activation products were created at negligible levels. All 
samples from In'adiation B exhibited no significant change in the light yield relative to the control 
samples. Their average dose and relative light yield is given in the figure, with the dose uncertainty 
covering the range of doses received by these samples (20 Gy to 90 Gy). The light yield uncertainty 
for all samples is the combination of the measurement standard deviation and an estimated relative 
uncertainty on the bias correction of 50%. A small loss of light yield was observed at approximately 
320 Gy for both WbLS and liquid scintillator samples. For doses of approximately 800 Gy, a clear 
degradation of the scintillation yield is evident in both samples. This damage level agrees with the 
online measurements within their uncertainty, suggesting that there was no measurable recovery of 
scintillator in the time between the irradiation and the measurement. The scintillation light yield 
was measured repeatedly for approximately one week for all samples following irradiation, and no 
recovery or further deterioration was observed. 
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Figure 5. The measured light yield in the liquid scintillator (red) and Water-based Liquid Scintillator (blue) 
following proton irradiation, relative to an unirradiated control. 

The optical absorbance spectra (figure ^ of both samples exhibit an increase with dose at 
wavelength regions relevant for the detection of scintillation photons. The reduction in the mea¬ 
sured light yield of this increased optical absorption can be quantified using an effective absorbance 
due to the radiation damage, oCgff'. 

fLY = e-^^^ff (3.1) 

where fu is the factor by which the light yield is reduced, and jc is the mean optical path length of 
the emitted photons within the scintillator. The effective absorbance is calculated by weighting the 
excess damage-induced absorbance (a(A), the difference between the irradiated and unirradiated 
traces in figure ^ with the normalised scintillator emission spectrum (S(A)) and photomultiplier 
quantum efficiency {QE{X))\ 

a,ff = J a{X)S{)i)QE{X)dX (3.2) 

The amount of the observed light yield change that is due to the damage-induced optical ab¬ 
sorbance is given in table approximating Jc to the radius of the sample vial. The additional 
optical absorbance induced in the samples is insufficient to fully account for the observed change 
in the scintillation light yield, which suggests that the damage is also evident as a reduction in the 
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scintillation photon production. This may be due to the radiation-induced formation of chemical 
quenching agents, or modification of the solvent/fiuorophore molecules involved in the scintilla¬ 
tion process um. We performed Fourier transform infrared spectroscopy to identify chemical 
changes in the samples. No significant changes were observed. Given that the doses in Refs [0, HO] 
were orders of magnitude larger than the maximum dose of this work, this is perhaps unsurprising. 


Table 2. The effect of damage-induced optical absorbance upon the measured light yield. 


Material 

Mean Dose 

Measured lighf 
yield loss 

Lighf yield loss expected 
from opfical absorbance 

Liquid Scinfillafor 

327 ± 10 Gy 

0.7 ± 0.4 % 

0.1 ±0.2% 

Liquid Scinfillafor 

797 ± 24 Gy 

1.7 ±0.3 % 

0.3 ± 0.2% 

Wafer-based Liquid 
Scinfillafor 

342 ± 10 Gy 

0.8 ± 0.5 % 

0.0 ± 0.2% 

Wafer-based Liquid 
Scinfillafor 

834 ± 25 Gy 

1.3 ±0.5 % 

0.7 ± 0.2% 


4. Discussion and Conclusion 

No damage to the scintillation light yield was evident for proton doses up to 90 Gy. Given that 
a typical patient dose in a proton therapy treatment fraction is 2 Gy, this suggests that a liquid 
scintillator-based proton therapy dose-verification system could service at least 45 patients with no 
change in performance. 

We will now consider the performance of such a QA instrument in more detail, given the ob¬ 
served changes at 800 Gy of proton dose. The typical annual number of patients served by a single 
proton therapy treatment room is approximately 300. Therefore, a QA instrument that validates the 
treatment plan of every patient may receive approximately 600 Gy in a year of operation, assuming 
that the instrument is used to validate a single 2 Gy treatment fraction for each patient. While our 
study has shown that the scintillator is damaged by this level of dose; an important consideration 
is that the entire volume of the scintillating phantom is not in'adiated during the QA measurement. 
Indeed, for a phantom volume of 30^ cm^, and a typical lateral target size of 9 cm^ and tai'get depth 
of 20 cm, the irradiated volume is approximately 1% of the total phantom volume. We assume a flat 
dose profile delivering uniform dose fo fhis volume, which is a conservative approximation. There¬ 
fore, for complefe mixing of fhe irradiafed scintillator wifh fhe larger volume, fhe lighf yield damage 
is dilufed by a facfor of 100 relafive fo our measuremenfs. In our measuremenf geomefry, 600 Gy 
of profon dose to a scinfillafor resulfs in an apparenf lighf yield loss of approximafely 1.3%. Half of 
fhis loss is assumed fo be due fo damage-induced opfical affenuafion increases, which corresponds 
fo Ueff = 0.0083 cm^^ Therefore, fhe irradiafed phantom volume would exhibif a scinfillafion 
lighf production reduction of 1.3%/(2 x 100) = 0.007%. For a mean optical pafh lengfh of 15 cm, 
fhe lighf absorpfion in fhe scinfillafor volume would decrease fhe defecfed number of photons by 
1 15*0.0083/100 _ 0 1 % Therefore, fhe amounf of lighf yield reduction due fo radiation damage 







in a QA phantom will be approximately 0.1% over a year of full-time clinical use. It is probable 
that this level of systematic variation of a dose verification instrument is acceptable for clinical 
operation given the magnitude of other uncertainties in the dose measurement flTTl] . 
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